The positron density in the intergalactic medium and the galactic 511 keV hne 
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The 511 keV electron-positron annihilation line, most recently characterized by the INTEGRAL/SPI 
experiment, is highly concentrated towards the Galactic centre. Its origin remains unknown despite 
decades of scrutiny. We propose a novel scenario in which known extragalactic positron sources such 
as radio jets of active galactic nuclei (AGN) fill the intergalactic medium with MeV e'^e~ pairs, 
which are then accreted into the Milky Way. We show that interpreting the diffuse cosmic radio 
background (CRB) as arising from radio sources with characteristics similar to the observed cores 
and radio lobes in powerful AGN jets suggests that the intergalactic positron-to-electron ratio could 
be as high as ~ 10""", although this can be decreased if the CRB is not all produced by pairs and 
if not all positrons escape to the intergalactic medium. Assuming an accretion rate of one solar 
mass per year of matter into the Milky Way, a positron-to-electron ratio of only ~ 10~^ is already 
enough to account for much of the 511 keV emission of the Galaxy. A simple spherical accretion 
model predicts an emission profile highly peaked in the central bulge, consistent with INTEGRAL 
observations. However, a realistic model of accretion with angular momentum would likely imply a 
more extended emission over the disk, with uncertainties depending on the magnetic field structure 
and turbulence in the galactic halo. 



The long-standing problem of the origin of the Galac- 
tic 511 keV line remains unsolved, four decades after its 
first detection. Observations by the INTEGRAL/SPI ex- 
periment imply the annihilation of at least ^ 2 x lO'*"^ 
positrons per second, mainly from a spherical bulge at 
the galactic centre subtending an angle of ~ 10° in the 
sky [1] [2] . The bulge-to-disk ratio of luminosities from 
e+e" annihilation is ~ 1.4 [3] and the ratio of the 511 
keV line to the lower energy continuum is consistent with 
100% annihilation via positronium formation (!]. 

Positrons are copiously produced in astrophysics. Stel- 
lar nucleosynthesis and supernova explosions produce ra- 
dioactive nuclides that decay via /?+ emission, resulting 
in MeV-scale positrons. High-energy particles, either 
photons or cosmic rays, can produce >> MeV positrons 
via pair production or pion decay. The rate of positron 
production in the Milky Way has been estimated for a 
number of candidate sources, but none reproduce the 
correct morphology (see [5j for a comprehensive review). 
Positrons produced by these sources can travel long dis- 
tances before annihilating, reducing the correlation be- 
tween the distribution of sources and the detection pat- 
tern, although this does not lead to any concentration 
of the signal in the bulge [B] . In light of the difficulties 
to explain the properties of the Galactic 511 keV emis- 
sion with astrophysical sources, annihilation or decay of 
dark matter has been invoked as an alternative source of 
positrons and as an explanation of the large contribution 
from the bulge {e.g. U^W)- 

This paper discusses a novel possibility. Positrons that 
escape from their source and their host galaxy in jets 
or winds may reach the intergalactic medium (IGM) to 
stay there indefinitely, as long as their cooling time re- 



mains longer than the age of the universe. Positrons of 
very high energy {E > 100(1 + z)-^/^ MeV at redshift 
z) are slowed by Compton scattering with the CMB, 
and positrons of very low energies (highly subrelativis- 
tic) cool by Coulomb scattering with thermal electrons. 
At intermediate energies, however, positrons do not cool 
over the age of the universe in the IGM. The turbulent 
magnetic structure of the escaping jets and winds can 
furthermore trap the pairs, forcing them to follow the 
motion of the plasma fluid as long as the magnetic field 
has enough small-scale structure to prevent particles from 
diffusing over large distances by moving along the field 
lines. Positrons fill the IGM, and if they are eventually 
accreted by gravity wells such as our own galaxy, they an- 
nihilate whenever they reach a medium with high enough 
density to cool by Coulomb scattering (e.g., [S])- 

To the best of our knowledge, the intergalactic positron 
abundance produced by astrophysical sources has not 
been previously considered. In this letter, we estimate 
the intergalactic positron density and show that positrons 
accreted by the Milky Way may substantially contribute 
to the Galactic 511 keV emission line. We first estimate 
the IGM positron-to-electron ratio based on cosmic radio 
background (CRB) observations, assuming that the ra- 
dio signal is produced in sources similar to the observed 
active galactic nuclei (AGN) jets and radio lobes. We 
then compute the expected 511 keV line luminosity and 
the emission profile in a spherical accretion model. We 
find that the annihilation in this spherical model happens 
mostly in the galactic bulge as required by observations. 
We conclude with a discussion of the expected changes 
of the emission profile in a more realistic accretion model 
as well as the potential of new radio observations of the 



small scale fluctuations to test our hypothesis. 

The CRB has been measured from 0.01 to 90 GHz 
by several collaborations including, most recently, the 
ARCADE-2 experiment [TU]. The observed brightness 
temperature spectrum, in excess of the constant cosmic 
microwave background (CMB) and corrected for Galactic 
emission, is well fitted by a power-law: 
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The shape of the spectrum is reproduced by syn- 
chrotron emission [TTJ [T^] from a population of elec- 
trons and positrons with a power-law distribution of 
energies E = jrUeC^, np{'j) dj ex j~p dj, with an in- 
dex p = 2.2. Extrapolations from luminosity functions 
of known synchrotron-emitting sources account only for 
about one sixth of the observed background intensity |llj . 
Underestimating the level of Galactic emission is a po- 
tential contaminant [T3]. However, the expected contri- 
bution is determined with tight errors, 5 (0.4) mK at 3.3 
(10) GHz, compared to the CRB brightness temperature 
of 54 (3) mK at the same frequencies. 

The brightness temperature in equation (II]) corre- 
sponds to an energy density per unit frequency, u^, of 
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= (1.0±0.1)xlO-^(^^) ■ eVcm-3 (2) 

Pairs with a Lorentz factor 7 radiating synchrotron 
in a magnetic field B, with Larmor frequency ujb = 
eB/{mec), dominate the emitted power at an emitted 
frequency v^ = Aj'^ujb and have a synchrotron cooling 
time tc = 9c/(47rea;^). Here, re is the classical elec- 
tron radius and A is a constant close to unity, which we 
use to match the integrated energy density in pairs to 
the photon luminosity at a given frequency (see equation 
4 below). We assume an approximately constant mag- 
netic field, and we define ie as the time over which pairs 
in a radio source radiate via synchrotron emission. This 
roughly corresponds to the lifetime of the source. During 
this time, pairs with Lorentz factor 7 will have radiated a 
fraction te/tc of their rest mass energy. This may be used 
to relate the average number density of radiating pairs 
with energy ^nieC^ at a mean source redshift (1 -|- Zr) 
to the present comoving energy density of emitted syn- 
chrotron photons: 
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where // is the fraction of contributing pairs that were 
able to escape to the IGM, and the observed frequency 
is v — I'e/il + Zr). Note that te/tc is much less than 
unity for low-7 pairs that radiate only a small fraction of 
their energy to the radio background, and can be larger 



than one for very high 7, when pairs need to be reacceler- 
ated many times within a source in order to maintain the 
emitted power-law spectrum. Now, using u^ ex i/^(p^^)/'^ 
and 7 — [z^(l -I- Zr)/{Au}B)]^^'^, and computing the to- 
tal number density of pairs Np — np{-f)jP/[{p — 1)7^~„] 
integrated for all 7 > 7min, we infer from equation ([3| 
that 
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This relationship tells us that the number density of pairs 
is equal to the energy density per unit log v in the ra- 
dio background divided by the electron rest-mass energy, 
times a number of dimensionless factors. Note that Np 
is independent of the frequency v at which we choose to 
evaluate it, because vu,, ex jy(3-p)/2^ rpj^g uncertainties 
are the magnetic field, the emission time ig, the mean 
source redshift Zr and the escape fraction //. The factor 
7„iin should be close to unity in the source rest-frame 
for a realistic mechanism to accelerate the synchrotron 
emitting particles. 

Using typical values for the lobes of radio-loud AGN 
of B ~ 10 fj,G and ig ~ 10^'^ yr (assuming that radio 
lobes expand at the typical intracluster sound speed), 
replacing p = 2.2 and using Zr ^ 1 and A = 1, and 
dividing by the mean electron density in the universe at 
present, fie = 2.3 x 10~''cm~^ [13], we obtain the IGM 
positron-to-electron ratio: 
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Thus, for typical parameters of a radio lobe, Np/fie — 
10~^//. For the core regions of AGN jets there is an ad- 
ditional Doppler factor of ^ = [r(l — I3cos0)]~^ in the 
observed luminosity due to relativistic beaming, where 
r is the jet bulk Lorentz factor and 9 is the angle be- 
tween the jet and the line of sight. Equation (Jsl is then 
modified with a factor 5 s^ . The typical value of the 
magnetic field in the core of AGN jets \s B ^ 3mG, as 
estimated from a sample of resolved sources [15] showing 
a characteristic spectrum of synchrotron radiation and 
self- absorption [TH]. For the same sources (5 ~ 8.5 and 
the emission time is the light-crossing time of the core re- 
gion, te ^ 25 yr. We obtain the IGM positron-to-electron 
ratio Np/fie ^ 10~^ //. Note that many radio jets may 
be produced in isolated galaxies without a massive intra- 
cluster medium surrounding them, so the relativistic jet 
may be directly expelled to the IGM without producing 
observable radio lobes. 

We conclude here that if most of the CRB were gen- 
erated by electron-positron pairs in sources similar to 
the AGN jets and radio lobes (with similar values of 
B^-^te), and if most of the positrons contained in these 
sources eventually escaped to the IGM, then there would 
be about ten positrons for every million electrons in the 



universe. This fraction may be much lower if the ra- 
dio background is mostly produced by electrons rather 
than pairs, or from sources with high values of B^'^t^. 
compared to AGN, or if most of the positrons do not es- 
cape. We note that these energetic particles would con- 
tribute to the matter pressure of the IGM. The ratio of 
the energetic particles pressure to the thermal IGM pres- 
sure is ~ NprUeC^ / {fiekBTi) ~ 5, for T/ = 10"' K and 
Np/Ue ~ 10~^, SO for the maximum value of the positron 
abundance the pressure would appreciably modify the 
IGM dynamics determining the properties of the Lya 
forest. 

Such a large e"*" fraction could lead to a visible annihi- 
lation signal as these positrons are accreted into galaxies 
with other intergalactic matter. In fact, if our galaxy is 
accreting matter at a rate M, the number of positron an- 
nihilations implied is lO**^-^ (iVp/rie) M/(M0 yr'^) s"^ 
For the Milky Way, we can expect a gas accretion rate 
of ~ lAf0yr~^ [T7], and reproducing the observed IN- 
TEGRAL 511 keV luminosity implying an annihilation 
rate of 2 x 10^"^ s^^ therefore requires Np/fif. ^ 10^^, 
which is safely below our estimated maximum and can 
be matched by assuming // ^ 10~^. 

Next, we model the distribution of the 511 kcV emis- 
sion from spherically symmetric accretion of plasma con- 
taining relativistic positrons. The dominant contribution 
to the positron cooling rate is due to Coulomb scattering, 
and is roughly proportional to the density of electrons in 
the interstellar medium (at high positron energies, the 
cooling rate does not depend much on whether the elec- 
trons are in ionized, atomic or molecular matter). A 
positron with an initial energy i^icM before infall will 
annihilate when it reaches a region of high enough den- 
sity for it to efficiently thermalize; this location depends 
on the matter distribution and the energy i^icM- 

We parameterize the matter distribution in the Milky 
Way as an axially symmetric distribution of matter, us- 
ing a model that compiles the results of several observa- 
tion studies described in detail in Ref. [TS] . The height- 
dependence of the Hi and H2 densities is Gaussian, with 
scale heights of 250 pc and 70 pc, respectively. The Hi 
density on the galactic plane at our location is ^ 1 atom 
cm~'^, whereas H2 is in a molecular ring around a galac- 
tocentric radius of i? = 5 kpc, with a peak density on 
the plane of ^ 2.5 atom cm^'^. The Hii distribution has 
two components, a central exponential disk with a scale 
height of Ikpc and central density of 0.025 atom cm~'^, 
and an annulus around R = 2 kpc representing Hii re- 
gions, with a scale height of 150 pc and peak density of 
0.2 atom cm~^. We use the energy loss-rates computed 
with the GALPROP package [Ull^, which also includes 
the subdominant effects of inverse-Compton scattering of 
CMB and interstellar light, bremsstrahlung, synchrotron 
radiation and ionization. Details of the exponential mag- 
netic field model and the interstellar radiation field used 
to compute these are also in Ref. [151 and references 



therein. 

The positrons are initially distributed in a power law 
with index of —2.2, following the inferred synchrotron- 
producing spectrum from the CRB, with lowest energy 
of 1 MeV. We have checked that decreasing the index 
to —2.5 (as in spectra observed in most AGNs) has lit- 
tle quantitative impact. Note also that positrons of very 
high energy should have cooled below ~ 10 MeV in the 
IGM by Compton scattering, although these are a mi- 
nority of the particles and it also has little impact on 
the computed annihilation distribution. The calculation 
of energy loss is stopped when positrons reach energies 
below 100 eV, at which point thermalization and annihi- 
lation should quickly follow. The flux of 511 keV photons 
can then be found from the annihilation rate per unit vol- 
ume dNf,+ /dVdt: 
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where x is the direction along the line of sight (l.o.s.), 
the prefactor of 2 accounts for the two outgoing pho- 
tons per annihilation and the 1/4 accounts for the fact 
positronium decay to two photons only occurs when para- 
positronium is formed, one quarter of the time. 

The results are then compared with the INTE- 
GRAL/SPI observations. Figure IT] shows a comparison 
between measured and predicted flux in the inner 16 
degrees of latitude (upper panel) and longitude (lower 
panel). These figures illustrate that the spherical ac- 
cretion model allows us to reproduce the large annihi- 
lation signal from the bulge. Using the parametrization 
of Ref. [2], we obtain a bulge-to-disk ratio of 2.6, leaving 
room for annihilations from Al^^ produced in supernovae 
and other /3"'"-producing isotopes distributed in the disk. 

Spherical accretion concentrates the infalling positrons 
close to the Galactic centre, producing a bright annihi- 
lation region in the bulge. A more realistic model would 
allow for the IGM gas to carry angular momentum and 
for a non-spherical galactic potential, so the gas could 
then fall over a much larger area around the disk. This 
would redistribute the luminosity over the disk, reducing 
the central luminosity and smearing the small-scale fea- 
tures. The positrons annihilate at the point where the 
medium is dense enough for them to lose sufficient en- 
ergy, and this may often occur in disk crossings rather 
than near the bulge. However, the details depend on the 
way the structure of the magnetic field in the Galaxy halo 
may modify the trajectories of accreting IGM plasma. 
Previous studies have considered a possible focusing ef- 
fect of a dipole magnetic field lines of the halo, although 
there has not been direct observation of such a dipole 
component [5j [21]. Nevertheless, if such a dipole field 
is present, it could contribute to an enlarged positron 
density in the bulge. This is clearly a complex problem 
that we cannot address in this letter. We note also that 
if the extragalactic positrons contribute substantially to 



the total annihilation signal but in a more extended fash- 
ion, the concentration to the bulge might also be caused 
by recent AGN activity from the central black hole in 
the Milky Way having delivered positrons to the bulge 
region. 
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FIG. 1: The binned INTEGRAL/SPI 511 keV data from ^31, 
integrated over Ihatitudo] < 8° (|Longitude|< 8°), is shown by 
the data points in blue in the upper (lower) panel. The outer 
four bins are 16° wide, and inner bins are 3.2° wide. Black 
solid line shows the predicted flux in the spherical mfall model, 
binned in the same way. The red dashed line shows the shape 
of the predicted flux, with its area matching the normalization 
of the histogram. 

The sources of the CRB are still unknown. It has re- 
cently been argued [22] that the RMS fluctuations of the 
CRB are more than an order of magnitude smaller than 
those of the infrared background. If cosmological, the 
sources for the radio background might be either at high 
redshift or spatially extended (on Mpc scales), so that 
the clustering amplitude is substantially lower. Alterna- 
tively, the radio sources might not be highly biased, for 
example if they arise from a population of low-luminosity 
AGN in low-mass halos, and this would reduce their clus- 
tering. The new generation of radio observations, mea- 
suring the small scale fluctuations of the background, 
would be able to test whether the cosmological popula- 
tion of AGN cores and lobes are sourcing the radio back- 
ground. We note that there might be positrons injected 
in the IGM at redshifts higher than ~ 6, which would 
Compton-cool to non-relativistic energies in the IGM be- 



fore they are accreted to a galaxy, so that at present they 
would annihilate at lower gas densities leading to a 511 
keV emission extending to higher galactic latitude than 
observed by INTEGRAL. 

In summary, we have presented a new possible source 
of positrons that can contribute to the 511 keV emis- 
sion from the Milky Way: the accretion of intergalactic 
material containing positrons that have been produced 
over the age of the universe in high-luminosity objects 
that can expel matter to the IGM, such as the jets of 
active galactic nuclei. Based on the radio background 
intensity, we have obtained their e contribution to the 
IGM using Eqs. Q and ^. If the CRB conies from 
these synchrotron-emitting sources, we estimate a maxi- 
mum density of 10^^ positrons per electron in the IGM. 
The positron fraction required to produce the observed 
INTEGRAL 511 keV luminosity is ~ 10-^ weU below 
our maximum estimation. The main feature of the IN- 
TEGRAL/SPI observation, namely the large density of 
positrons in the Galactic bulge, can be reproduced by a 
simple spherical model of accretion, although a predic- 
tion of the exact line emission morphology will require 
more realistic modelling of IGM gas accretion and of the 
magnetic field structure in the Galactic halo. 
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